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Mapping tlie expansion of the universe gives clues to the underlying 
physics causing the recently discovered acceleration of the expansion, and 
enables discrimination among cosmological models. We examine the utility 
of measuring the rate of expansion, H{z), at various epochs, both alone and 
in combination with distance measurements. Due to parameter degeneracies, 
it proves most useful as a complement to precision distance-redshift data. Us- 
ing the baryon oscillations in the matter power spectrum as a standard rod 
allows determination of H(z) / (ílmh'^Y^'^ free of most major systematics, and 
thus provides a window on dark energy properties. We discuss the addition 
of this data from a next generation galaxy redshift survey such as KAOS to 
precision distance information from a next generation supernova survey such 
as SNAP. This can provide useful crosschecks as well as lead to improvement 
on estimation of a time variation in the dark energy equation of state by 
factors ranging from 15-50%. 



1 Introduction 

We now have strong evidence that the expansion of the universe is accel- 
erating, from the original method of Type la supernova distance-redshift 
measurements |2¡ and concordant observations of the cosmic microwave 
background (CMB) power spectrum and of large scale structure The 
nature of the dark energy responsible for the acceleration will have profound 
implications for cosmology, particle physics, and fundamental physics. Map- 
ping the expansion history of the universe offers a way to gain insights into 
the dark energy and the fate of the universe, for example by characterizing 
the equation of state behavior that is directly related to properties of the 
scalar field potential. 

As discussed in Linder one would like to carry out this mapping with 
not only precision measurements of the distance-redshift relation, but ide- 
ally with data on differential distances corresponding to the change between 
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neighboring redshift epochs. The former, notably from the Type la supernova 
method, have proved adept at constraining the energy density and equation 
of state of the dark energy, with great improvements expected in the next 
decade. But these involve an integration over the expansion rate behavior 
H{z), which itself involves a redshift integral over the equation of state w{z). 
Probes more closely related to the differential distance might give H{z) more 
directly. 

However the integral nature of the distance-redshift relation also provides 
the power to break degeneracies between cosmological parameters, which is 
an equally important aspect. So [Sj found that the Alcock-Paczyhski effect 
of the cosmic shear distortion - due to the source distances radial and trans- 
verse to the hne of sight being measured at different epochs - did not in 
fact automatically give more stringent estimations of the dark energy prop- 
erties, despite involving a bare factor H{z). The cosmic shear (not to be 
confused with the local, weak lensing shear) is related to the ratio of the 
differential distance over some redshift interval to the integrated distance to 
the source. So it is interesting to consider whether the situation changes if we 
can independently measure the two quantities, basically finding the Hubble 
parameter H{z) separately. 

In Section|21we investigate the use of H{z) for constraining the cosmolog- 
ical model. But in SectionlSlwe find that the most promising technique - the 
baryon oscillation method - actually measures a slightly different quantity. 
We then examine the use of the radial and transverse distances provided by 
precision next generation galaxy redshift survey observations of the linear 
matter power spectrum, separately and together. In Section IU we show that 
the fuU power of the method comes from adding the information to a deep 
distance survey such as from accurate observations of Type la supernovae 
(e.g. SNAP). We summarize our conclusions and the need for future work in 
SectionEl 

2 Using H{z) Information 

In this section we consider a data set giving the Hubble parameter H{z) 
at some redshifts z, with a certain fractional precision. This is a purely 
theoretical investigation as we do not specify how the measurements are 
made. Indeed, as mentioned in the introduction, the cosmic shear method 
only gives the product of H{z) with the distance corresponding to the redshift 
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and as we will see in ^the baryon oscillation method also provides a ratio 
involving H{z). So this is meant as a thought experiment. 

Similarly, it is obvious that knowledge of H{z) over the entire redshift 
range from the observer at 2; = out to some depth is overly optimistic and 
would supersede any distance measurements in that range. So we consider 
data at only a few redshifts in a narrow range and ask what cosmological 
information this can provide and what value it adds to a more reahstic set 
of distance measurements. Recall that the comoving distance r or conformal 
time 1] is related to H{z) in a flat universe by 



and the angular diameter distance da = {l + z)~^r and luminosity distance 
di = {1 + z) r. The differential distance along the hne of sight (radially) is 
simply dr\\ = drj = dz / H and transversely is dr± = rjO, where 6 is the angle 
subtended. 

Through the Friedmann equations, the expansion rate H{z) is related to 
the cosmological components by 



where Hq is the Hubble constant, the present value H{z = 0), Qm is the 
dimensionless matter density, and w{z) is the equation of state of the dark 
energy. We can examine the impact of measurements of H{z) on determina- 
tions of the cosmological parameters through the sensitivities dH/dílm etc, 
achieving formal constraints through the Fisher matrix method |ü]. 

The sensitivities are shown in Figure^ with the parametrization w{z) = 
Wq + Waz/ {1 + z) of Linder |7| that allows robust treatment of the equation 
of state to redshifts greater than one. However, the sensitivities are not 
the whole story: degeneracies between the parameters play a major role in 
determining them. For example, while a 1% measurement of H{z) at 2; = 3, 
say, would apparently constrain Qrn to 0.06, wq to 0.14, and Wa to 0.3, this 
holds only upon flxing all parameters but one. In fact a measurement at a 
single redshift only gives an inflnite elhpsoid in the joint three dimensional 
parameter space. Even over a redshift range, such as H{z) to 1% at z = 2.8, 
3, 3.2, the uncertainties are uselessly large: a^flm) = 0.87, cr{wo) = 76, 
<y{'Wa) = 207. But because the elhpsoid is fairly narrow, and the degeneracy 
direction is different than for distance measurements, the combination of 
H{z) information with distance information can be valuable. 




(1) 



{H/Ho)' = nm{l + z)' + (1 - fi^)e3r'^''^(i+^')[i+"'(^')l. 



(2) 
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Figure 1: The sensitivity of the expansion rate H{z) to the cosmological 
parameters is plotted as a function of redshift. The larger the derivative 
at a particular redshift, the more constraining the observations are there, 
but the curves by themselves do not account for degeneracies between the 
parameters. 



For example, adding the estimation of H{z) aX z = 2.8, 3, 3.2 to a simu- 
lation of the data expected from the Supernova/Acceleration Probe (SNAP; 
|H1) survey out to z = 1.7 allows parameter determination to <7{Qrn) = 0.0082, 
<j{wq) = 0.078, cr{wa) = 0.45. This represents a factor 3.5 improvement in 
constraining í]^, 2% in Wq, and 23% in the measure of the time variation 
w' ~ Wa/2, relative to the canonical SNAP results. So as expected there 
is clearly value in obtaining measurements of H{z) (though we have not es- 
tabhshed how such would be carried out) - though only in complementarity 
with a distance probe. 
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Indeed one can show that measurements of H{z) at redshifts z > 1 basi- 
cally act hke information about the matter density Qm- If one eschewed any 
H{z) data but added a prior an^ = 0.0082 to the SNAP data then one would 
roughly recover the previous parameter estimations. This is not surprising 
since at 2 > 1 one is increasingly in the matter dominated, deceleration epoch 
and the expansion rate therefore best measures the matter density, not the 
dark energy properties. So an integral measure such as the distance-redshift 
relation actually has an advantage in probing the dark energy equation of 
state, despite this quantity entering the distance through a double integraL 

We emphasize this important point further by considering two elabora- 
tions. If we spread the redshift range of the H{z) measurements, to z = 2.5, 
3, 3.5 and add simulated information from the future Planck cosmic mi- 
crowave background survey |Ü], then the dark energy constraints are still 
weak: a{ílm) = 0.039, a{wo) = 1.6, a{wa) = 5.4. Again, the CMB has 
limited sensitivity to the dark energy equation of state and little comple- 
mentarity with the H{z) measurement. If we now add the SNAP data, the 
estimations improve to 0.0056, 0.070, and 0.34 respectively, but httle of this 
is due to H{z) since the CMB complementarity is much stronger. The part 
of the improvement due to H{z) is mostly restricted to flm (since that is 
what H{z = 3) best probes) and somewhat wq (due to its degeneracy with 
fim); adding H{z) tightens estimation of flm by 44%, wq by 11%, but w^ by 
only 4%. 

For determination of H{z) near z = 1, the situation is only slightly better. 
It no longer acts as predominantly a matter density prior, but again H{z) by 
itself cannot constrain the dark energy parameters, even with observations 
over a range z = 0.5 — 1.5. Furthermore, it has less complementarity with 
SNAP data and improves w^ constraints by only 6%. But conversely it gains 
in complementarity with the CMB information, and improves SNAP+CMB 
parameter determinations of Qm, wq, Wa by 12%, 23%, 16% respectively. 
These various cases are iUustrated in Fig. |21 

So the distance data plays a central role in determining dark energy prop- 
erties and H{z) measurements only a subsidiary, complementary one. But 
in any case we note that we have not identified any cosmological probe that 
provides these Hubble parameter measurements, let alone at 1% accuracy. 
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Figure 2: Joint probability contours (68% confidence level) for the time vari- 
ation Wa and the present value of the dark energy equation of state show 
that measurement of ií{z) is not as powerful as the CMB in complementing 
supernova data. 

3 Baryon Oscillations 

One cosmological probe that shows promise in obtaining differential distance 
measurements is the use of the imprint of the primordial baryon-photon 
acoustic osciUations in the matter power spectrum. These are analogous 
to the oscillations appearing in the CMB temperature power spectrum, but 
are much smaller, appearing as wiggles superposed on the larger dark matter 
component (see Pül E] for a comprehensive discussion). The wiggle wave- 
length can be used as a standard ruler, since the intrinsic scale is known from 
well understood physics of the matter-radiation decouphng epoch. Then the 
angular or redshift space scale can be measured through a wide field red- 
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shift survey (though beyond the current state of the art) and the comparison 
probes the cosmological model. 

By observing at redshifts z > 1 some of the wiggles appear in the hn- 
ear density regime of the power spectrum, and by using only the locations 
and not the amphtudes of the osciUations one does not require problematic 
models of structure formation and evolution. This method then has several 
positive aspects: simple, hnear physics free from astrophysical uncertainties, 
direct relation of observations to cosmological quantities, and sensitivity to 
a snapshot of the expansion rate, H{z). For further discussion of the details 
and possible implementation of this probe see |T2l 1131 • 

However, the baryon oscillations do not provide a pure measure of H{z). 
Rather, the physics involves the ratio of the "standard rod" size to the ob- 
served osciUation scale (generally in Fourier wavenumber, fc-space). So the 
central quantity is 

K ^ — = -dz— = (3) 

kobs s dz H{z)s' 

where is the predicted acoustic osciUation scale, proportional to the inverse 
of the sound horizon s, and kobs is the observed scale, proportional to the 
inverse of the standard rod length d-q. The sound horizon is given by 



dz{cs/H) (4) 

= (a^/i')-^/' I dacs [a + a,, + (1 - fi-i)a"e-=^/'^i'^"[i+-(")]] (5) 

where a = {1 + z)^^ is the scale factor of the universe, Cg is the sound speed 
in the baryon fluid, z¿f.i. is the redshift of decouphng, and is the scale 
factor at matter-radiation equahty. Note that while a cosmological constant 
{w = —1) would cause the last term in the brackets to have a neghgible 
contribution to the integrand, some forms of dynamical dark matter could 
have non-neghgible influence at these redshifts (see |15^ for a discussion of 
early quintessence). The sound speed for adiabatic perturbations in the 
baryons is 

Since the baryon density pb ~ VLbh'^ is weU determined by current CMB 
measurements, and wiU be further improved by Planck data, and the photon 
density ~ is also accurately known, then we can regard Cg as flxed. 
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From the form of Eqs. Q and (@)), we see that H{z) enters in both 
numerator and denominator, as itself and as an integrand. This is the same 
form as for the cosmic shear probe [3]. So as pointed out there, the value of 
the Hubble constant Hq or h does not enter the problem and therefore does 
not require marginahzation. This is a definite advantage. Furthermore, one 
can divide numerator and denominator by (ílmh'^Y^'^. At a casual glance, 
one might think that all dependence on this quantity is then removed. But 
in fact, the approximation that s ~ is not a good one, as pointed 

out in ^ü]. There, a closer approximation for a fiat, cosmological constant 
universe was found to be s ~ (^2^^^)"° '^, while a more precise analysis |TT¡ is 
equivalent to s ~ (Í7m^^)~°'^^- The additional factors come from the presence 
of üeg and to a much lesser extent adec- However, since the dependence arises 
from the sound horizon, it is the same for all the redshifts at which the 
measurements of kobs are carried out by the redshift survey. This, combined 
with the precision to which Planck wiU determine ílmh'^, means that its 
uncertainty couples very weakly to the other parameters (this was exphcitly 
tested), and we wiU neglect it. However, one stiU must incorporate the 
uncertainty in ílrn in order to obtain reahstic parameter error estimations. 

Therefore, the baryon oscillation method can essentially provide measure- 
ments of two cosmological variables according to Eq. Q: H{z) = H{z) / {ytmh^Y^'^ 
and f]{z) = r{z) (íí^/i^)!/^, These come respectively from the wavenumbers 
along {di]) and transverse (77) to the hne of sight. This distinction from a 
plain H{z) as treated in ^is important for the parameter degeneracies and 
complementarity with other methods. 

Figure El shows the sensitivities for these two quantities. As expected, 
for wq and Wa the derivatives are the same as for H{z) and r{z). However, 
the degeneracy relations between the parameters have now changed, and so 
the strength of the estimations have as well. Again we find that the probe 
in isolation cannot effectively constrain the cosmological model - even the 
matter density since most of its dependence has been removed in the ratio. 
Even in combination with CMB data it has httle leverage. 

However the situation changes significantly for a fiducial model that has 
time variation in the equation of state. For a supergravity inspired model 
[T7¡ that is well fit by wq = —0.82, Wa = 0.58, the oscillation data offers 
definite sensitivity to the time variation. Now a 2% (1%) measurement of K 
in both its radial and transverse aspects, in combination with Planck data, 
allows estimation of w^ to 0.29 (0.20) for measurements a.t z = {0.5, 1, 1.5} 
and 0.62 (0.47) for z = {2.5,3,3.5}. However, the estimations of Wq remain 
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poor: 0.16 (0.09) and 0.36 (0.27) respectively. Furthermore, since only 1-2 
wiggles are in the hnear regime at 2; ~ 1, the observations are unhkely to 
achieve better than 2% precision there, while we see that even 1% precision 
at ^ ~ 3 gives less impressive results. And of course we have no guarantee 
that the true cosmological model wiU have a strong time variation in the 
dark energy equation of state. 




Figure 3: As Fig. ^ but for the baryon osciUation scale. Baryon osciUation 
measurements provide two quantities, corresponding to wave modes along 
[dr]] and transverse [1]] to the hne of sight. The flrn curves do not have zero 
sensitivity at z = because of the division by {ílmh'^Y^'^. 
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4 Baryon Oscillations Plus Supernovae 



As with the H{z) analysis in ^ the baryon oscillation method cannot stand 
alone as a robust cosmological probe. In this section we consider it in com- 
plement with a SNAP supernova distance survey. As expected, we find that 
it does not behave in the same manner as H{z), as effectively a prior on 
the matter density. In fact, unlike H{z), within a cosmological constant 
model of dark energy the complementarity with precision distance measure- 
ments is now substantial, providing good constraints. Adding oscillation 
information acts even slightly more strongly than adding CMB information, 
relative to supernovae. When the baryon oscillation information is added 
to SNAP-fCMB, further modest improvements are seen - around 14% in 
both wq and Wa for 2% measurement of the osciUation scale and 30% for 1% 
measurement. This is fairly insensitive to the exact redshift distribution of 
the matter power spectrum measurements, i.e. for redshifts near 1 or 3, or a 
spread z = {0.5, 1, 1.5} vs. {0.8, 1, 1.2}. 

For the SUGRA model, the improvement is stronger. Baryon oscillations 
and CMB have increasing complementarity to each other and to supernovae. 
Now upon adding the oscillation probe to SNAP-I-CMB the estimation of Wa 
sharpens by 46-37% (54-60%) for 2% (1%) precision, depending on whether 
the measurements are near z = 1 or 3. Furthermore, (t{wo) reduces by 51- 
39% (59-57%). So this offers hope that the baryon osciUation method can 
provide important complementarity useful in uncovering the nature of the 
dark energy. The error contours for the cosmological constant and SUGRA 
cases are shown in Figs. |3] and |S1 

Because of the insensitivity of the baryon oscillation results to the exact 
redshift range, so long as 2; > 1, one can choose the survey characteristics 
based on observational considerations. As mentioned, at redshifts z < 1.5 
one is likely to detect only 1 or 2 baryon wiggles, making it more difficult 
to precisely determine the osciUation scale kobs- But for z > 1.5 the linear 
regime quickly extends in fc-space, due to behavior of structure formation in 
a universe recently dominated by dark energy (see Fig. 1 of ^^), providing 
3-4 detectable osciUations. The redshift ranges most advantageous for ob- 
servations are often identified as 2; = 0.5 — 1.3 and z = 2.5 — 3.5 |T3] due to 
easy selection by 4000A break and Lya features in the galaxies used in the 
survey. Estimates of number of galaxies required and sky coverage are given 
in 01131. 

Such a redshift survey could be accomplished by large telescopes on the 
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Figure 4: As Fig. El for baryon osciUation measurements and a cosmological 
constant model. The baryon oscillation data is slightly stronger than the 
CMB data but not very complementary with it. 



ground within a decade. One possibihty is the KAOS project [TH]: the Kilo- 
Aperture Optical Spectrograph proposed as a front end for the Gemini South 
8 meter telescope. This would have multiplexing capabihty from some 4000 
fibers for simultaneous measurement of galaxy redshifts. With a 1.5 square 
degree field of view and coverage of some 400 square degrees of sky KAOS 
could measure precise redshifts for lO^ galaxies. This could provide estimates 
of the wiggle scale at the 2% precision level [Tü] . 

Another intriguing idea is to use wide field observations from space. This 
would have the advantage of not being restricted to the z ^ 1 and 3 ranges 
just mentioned, which were hmited by the Earth's atmosphere. Indeed, from 
a theoretical point of view, a redshift range z = 1.5 — 2 is essentiaUy as 
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Figure 5: As Fig. I^but for a supergravity inspired dark energy model. Now 
baryon osciUation and CMB data are more complementary with each other 
as well as the supernova data, so the improvement with all three is more 
pronounced. 



powerful as 2; 3 in terms of number of oscillations mapped, a definite 
advantage over z ^ 1, and yet requires less spectroscopic exposure time than 
the deeper survey. Calculations show that the parameter estimation for a 
given precision is as tight as the lower or higher redshift ranges. 

While there is no planned massively multiplexing spectrograph for space, 
one interesting possibihty is populating spare regions of the SNAP focal 
plane with grisms capable of low resolution spectroscopy. Also, photometric 
redshifts can be generated with SNAP's nine filters. There is no problem 
achieving the number or area statistics as the proposed SNAP wide field 
survey (mostly focused on weak gravitational lensing) wiU find lO^ galaxies 
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over 300 square degrees. However with a 2 meter telescope and spectral 
resolutions of order 100 or less, this clearly is not capable of carrying out 
all the science that an 8 meter ground based telescope with high resolution 
spectrograph could. StiU, while this would not provide the same precision 
mapping of the 3D matter power spectrum, it might give decent quahty 
information on the 2D, projected spectrum, roughly corresponding to the 
transverse wavenumber modes in Eq. Q. Photometric or low resolution 
spectroscopic redshifts would additionally give a smeared representation of 
the radial dimension. Detailed analysis of the baryon osciUation method 
with SNAP is left for future work; here we simply investigate the parameter 
constraints from the transverse and radial modes separately. 

One expects that the radial mode, dfj, which involves a bare factor H{z), 
should provide better hmits, while the transverse mode, f¡, acts basically 
hke a distance-redshift measurement though with the important degeneracy 
differences previously mentioned. Indeed the sensitivities plotted in Fig. El 
bear this out (though the degeneracy relations are not there apparent; also 
note that the radial mode has low sensitivity at redshifts that are well into 
the matter dominated epoch). For example, denoting the full baryon oscil- 
lation information as BO, the radial only as B0||, and the transverse only 
as BO_L, one finds that 2% precision gives <j{Qrn) = 0.0057, cr(wo) = 0.069, 
a{wa) = 0.30 for SN+CMB+BO, (0.0073, 0.073, 0.31) for SN+CMB+BO||, 
and (0.0065, 0.075, 0.35) for SN+CMB+BO^. (Though presumably the pre- 
cision in a fuU 3D survey would be better than in a 2D plus low resolution 
radial survey.) In the last case there is essentially no improvement over the 
SN+CMB case without any osciUation information. Various cases are iUus- 
trated in Fig. |ül 

Again, the baryon oscillation method is more useful in the presence of 
a time varying equation of state dark energy. Moreover, due to the long 
baseline entering the integrated distance, the BO^ information is actually 
more valuable than the radial, if at the same precision. The results for the 
SUGRA model are iUustrated in Fig. [7| We see that 2% precision data 
on the transverse modes adds complementarity to the supernova and CMB 
information, improving the estimates of Qrn by 56%, Wq by 46%, and w^ 
by 42%. This represents the vast majority of the impact of the full baryon 
osciUation data discussed at the beginning of this section. 
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Figure 6: As Fig. 0] but separating the contributions from the baryon os- 
cillation wave modes transverse and parallel to the hne of sight. Radial 
information, requiring accurate redshifts from a spectroscopic survey such as 
KAOS, is more usefuL 

5 Conclusion 

Differential distance measurements, providing a snapshot of the expansion 
rate H{z), have long seemed attractive theoretically as ways to probe the 
nature of dark energy. But they have also appeared difñcuh to implement 
observationally. The cosmic shear, or Alcock-Paczyñski, probe (not to be 
confused with promising weak lensing shear measurements) involves a ra- 
tio of differential to integrated distances, or the product H{z)r{z), and |3] 
showed that it could act only in a minor, complementary role to precision dis- 
tance observations. The growth of structure in the hnear regime also might 
be thought sensitive to H{z) but at redshifts z > 2 this essentially probes 
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Figure 7: As Fig. El but for a supergravity inspired dark energy model. Now 
the transverse modes, detectable through a 2D survey, are more influentiaL 



Í2m not the dark energy; however, through other factors the hnear growth 
stiU retains some sensitivity to the equation of state and its time variation 
¡21¡ . Nonhnear structure formation can involve H{z) as a separate factor 
through the differential volume element in cluster or galaxy halo counts, but 
this is entangled in systematic uncertainties from nonhnear physics and ob- 
servational selection effects fÜll2D|- On large scales this may be amehorated, 
but the needed numerical simulations of large scale structure incorporating 
a time varying equation of state are just now being carried out JTT . 

In this paper we have pushed these observational difñculties into the 
background and considered the use of H{z) regardless. Our conclusion is 
that it is not a panacea and only offers aid through complementarity with 
a deep, precision distance survey such as SNAP; then it contributes mild 
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improvement to the cosmological parameter constraints. This is basically due 
to H{z) acting at higher redshifts as a determinant of the matter density, 
not a direct probe of dark energy properties. At redshifts z < 1 it has 
somewhat more leverage, but requires precision on the 1% level for significant 
improvement. 

The baryon osciUation method of using wiggles in the matter power spec- 
trum as a standard ruler determines a shghtly different measure of the ex- 
pansion rate H{z). This probe is sufñciently promising, though again only in 
complementarity with a supernova distance survey, that it should be pursued 
further. Depending on the nature of the dark energy, incorporation of oscil- 
lation measurements can oífer significant improvements on estimation of the 
time variation of the equation of state. One of the most striking aspects is 
its cleanness, based on simple, weU understood physics and with no apparent 
major systematic uncertainties. Note that in all the analyses presented here 
of different cosmological probes, only the SNAP data has included systematic 
uncertainties - the H{z) and baryon osciUation precisions have been taken as 
purely statisticah For all known and other proposed probes this is certainly 
overly optimistic. Whether systematics enter at the 1-2% level in the baryon 
oscillation method, from, say, residual nonlinearities or mass vs. light bias, 
needs further investigation. 

Two interesting concepts for baryon osciUation observations are the KAOS 
project on the ground, and a spectroscopicaUy less precise but reasonably 
straightforward space implementation with grisms or photometric redshifts 
from SNAP. We have seen that the optimal redshift range is not strongly 
determined by the parameter sensitivity, and so wiU be driven by trade offs 
in observation strategy. Both projects deserve further investigation, though 
it is intriguing to imagine that SNAP could represent a cosmology super- 
probe - incorporating the supernova distance, weak lensing, some part of the 
baryon osciUation, and possibly even cluster count methods of cosmological 
parameter determination. But even if SNAP is rather promising for reveal- 
ing the nature of dark energy, our understanding and confidence wiU stiU be 
strengthened by multiple, complementary and crosschecking next generation 
surveys. 
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